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I.-  BACKGROUND;  REVIEW  OF  BASIC  CONCEPTS. 


This  fin*!  report  illustrates  the  resalts  of  a  three-years  research  program  on  Quan¬ 
tum  Chaos. 

After  some  ten  years  of  investigations  on  the  quantum  dynamics  of  classically 
chaotic  systems,  the  basic  question  whether  in  quantum  mechanics  anything  survives 
of  the  impressive  manif  stations  of  classical  chaos  hasn’t  yet  been  given  a  dear  answer. 
If  a  general  indication  is  to  be  drawn  from  the  analytical,  numerical  and  experimental 
results  up  to  now  obtained,  this  is  that  classical  chaos  is  suppressed  or  at  least  strongly 
inhibited  by  quantization;  as  a  consequence,  "quantum  chaos”  -  a  widespread  denomi¬ 
nation  for  this  research  area  -  is  still  sometimes  consider ered  a  questionable  concept.  In 
spite  of  this  negative  remark,  the  search  for  quantum  manifestations  of  dassical  chaos 
has  led  to  substantial  advance  in  several  areas  of  microphysics:  semiclassical  methods 
Lave  been  substantially  refined  and  new  phenomenological  aspects  .f  the  dynamics  of 
small  quant  tun  systems  have  been  identified. 

1.1-Chaos  in  microphysics? 

A  major  impulse  to  the  investigations  reported  below  was  given  by  the  analysis  of  some 
specific  models  used  in  various  areas  of  microphysics.  From  the  standpoint  of  classical 
mechanics,  these  models  exhibit  the  typical  feature:  of  nonlinear  Hamiltonian  systems 
and  under  appropriate  conditions  they  enter  a  chaotic  regime.  The  question  is,  how  is 
the  classical  chaos  mirrored  in  the  quantum  dynamics  of  these  models? 

This  is  a  very  important  question  for  the  present  day  physics,  both  from  a  the¬ 
oretical  and  an  applicative  viewpoint.  On  theoretical  grounds,  it  involves  the  nature 
of  the  semiclassical  approximation  when  the  classical  motion  is  chaotic,  i.e,  in  a  situa¬ 
tion  in  which  the  boundary  between  quantum  and  classical  mechanics  is  still  vague  and 
undefined  so  many  years  after  the  birth  of  quantum  mechanics. 


On  more  physical  grounds,  the  onset  of  classical  chaos  is  often  accompanied  by 


strong  instabilities  that  lead  to  experimentally  observable  and  often  dramatic  behavior. 
The  possibility  that  this  sort  of  phenomena  can  be  observed  in  microphyrics,  i.e.  in 
the  domain  of  validity  of  quantum  mechanics,  opens  a  vaste  new  field  to  experimental 
research.  A  widely  investigated  example  (  to  be  discussed  in  more  detail  below  ),  to 
which  the  work  reported  here  has  given  essential  contributions,  is  that  of  a  atom  in 
an  external  radiation  field.  This  involves  the  study  of  the  dynamics  of  an  electron  in 
the  combined  Coulomb  plus  radiation  field,  and  a  simple  model  for  that  is  a  nonlinear 
oscillator  subject  to  an  external  perturbation  periodic  .n  time.  Such  a  classical  model 
turns  out  to  exhibit  a  transition  to  chaos  as  soon  as  the  perturbation  strength  reaches 
a  critical  value.  Were  the  atom  a  classical  object,  in  such  conditions  it  would  be  rapidly 
ionized  (chaotic  ionization).  The  atom  however  is  quantum,  so  the  important  question 
arises:  does  quantum  mechanics  allow  for  any  such  ionization  mechanism? 

According  to  experimental  and  theoretical  results,  it  does  -  under  certain  conditions 
at  least.  Nevertheless,  the  modifications  that  quantum  mechanics  imposes  over  the 
classical  chaotic  picture  are  extremely  important  and  their  investigation  has  opened  a 
new  field  to  both  theoretical  and  experimental  analysis. 

1.2-Dynamical  localization 

The  concept  of  dynamical  localization  dominates  the  present-day  understanding  of  quan¬ 
tum  chaos  in  periodically  perturbed  systems.  A  classical  nonlinear  oscillator  subject  to 
an  external  perturbation  periodic  in  time  can  enter  a  regime  of  chaotic  motion  marked 
by  diffusive  growth  of  its  energy.  The  motion  of  such  an  oscillator  becomes  very  erratic 
and  looks  very  much  like  a  random  walk;  in  that  case,  only  a  statistical  description 
proves  viable  (  though  the  system  is,  in  principle,  perfectly  deterministic  and  has  just 
a  very  few  freedoms  ).  This  kind  of  a  motion  is  called  chaotic  diffusion. 

Chaotic  diffusion  is  a  physically  impoilant  aspect  of  classical  chaos.  Chaotic  in¬ 
stability  leads  to  this  type  of  response,  e.g  ,  in  nonlinear  systems  subjected  to  periodic 
perturbations.  This  type  of  instability  is  well  known  in  classical  dynamics  and  is  often 
responsible  for  undesirable  effects  in  macroscopic  systems. 


3 


The  stability  of  atoms  or  molecules  in  external  fields  clearly  depends  on  whether 
this  phenomenon  survives  at  least  in  part  at  the  quantum  level.  It  turns  out  that 
quantization  has  a  strong  effect  on  this  diffusion.  As  a  matter  of  fact,  the  study  of  a 
simple  albeit  abstract  model  (  the  so-called  Kicked  Rotator  )  yielded  the  remarkable 
result,  that  quantization  strongly  limits  and  even  suppresses  this  type  of  diffusion.  In 
other  words,  a  classical  system  would  continue  to  absorb  energy  from  the  external 
perturbing  field,  but  its  quantum  version  is  much  more  stable  and  reaches  after  a  while 
a  sort  of  steady-state  in  which  no  energy  is  further  absorbed  (  in  the  average  ). 

The  mechanism  that  produces  this  quantum  quenching  of  the  classical  diffusive 
behaviour  is  basically  quantum  interference. 

In  the  chaotic  regime,  the  classical  oscillator  would  display  a  host  of  different  orbits 
that  behave  quite  differently  in  spite  of  the  closeness  of  their  starting  points  ( initial  con¬ 
ditions);  in  the  quantum  case,  all  these  wildly  different  dynamical  possibilities  interfere, 
thus  leading  the  process  of  excitation  to  a  halt. 

This  phenomenon  has  been  recognized  to  be  analogous  to  the  Anderson  localization 
of  Solid  State  Physics.  1 

At  very  low  temperatures,  the  resistance  offered  by  a  cristalline  solid  to  the  mo- 

I 

tion  of  electrons  is  essentially  due  to  imperfections  in  the  crystal  structure,  that  spoil 
its  per  fee,  periodicity.  The  Anderson  model  is  a  model  fo*  the  study  of  the  residual 
resistance,  in  which  the  electron  jumps  from  one  site  to  another  in  a  disordered  lat¬ 
tice  -  i.e.,  a  lattice  that  exhibits  random  deviations  from  the  perfect  crystal  periodicity. 
Classically  a  random  walk  would  be  expected,  but  quantum  mechanically  it  turns  out 
that  the  diffusion  of  the  particle  through  the  disordered  lattice  is  stopped  or  at  least 
inhibited  by  interference  (  depending  on  the  spatial  dimension  of  the  lattice  ).  Due  to 
this  effect,  the  residual  conductance  of  a  linearly  shaped  sample  of  a  cristalline  solid 
decreases  exponentially  with  the  length  of  the  sample.  Ihe  rate  of  exponential  decay 
7  across  a  given  sample  depends  on  the  specific  sample  ;  in  the  Anderson  and  related 
models,  the  lattice  imperfections  are  described  by  a  random  potential,  so  that  the  actual 
rate  of  decay  depends  on  the  particular  realization  of  that  random  potential.  However, 
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in  the  limit  of  an  infinitely  long  sample  7  turns  out  to  be  realization-independent  and 
its  inverse  7-1  is  known  as  the  localization  length. 

Realizing  that  the  effect  of  quantization  on  classical  chaotic  diffusion  and  the  An¬ 
derson  localization  are  just  different  manifestations  of  a  common  quantum  dynamical 
effect  was  the  starting  point  of  a  theory  known  as  dynamical  localization  theory.  At 
the  heart  of  this  theory  lies  the  assumption  that  this  effect  is  a  generic  occurrence  in 
the  quantum  dynamics  of  systems  that  classically  would  display  a  chaotic  diffusion,  and 
that  it  can  be  efficiently  analyzed  by  means  of  the  same  methods  used  in  the  theory  of 
Anderson  localization.  This  means  that  the  dynamics  should  be  analyzed  on  statistical 
grounds.  In  other  words,  the  theory  should  give  up  a  complete  description  of  all  the 
details  of  the  quantum  dynamics;  only  certain  important  average  quantities  like  the 
localization  lengths  should  be  taken  into  account. 

The  link  of  the  dynamical  localization  effect  with  the  Anderson  localization  can  also 
be  qualitatively  understood  as  follows.  A  quantum  system  subjected  to  some  external 
monochromatic  field  will  perform  a  chain. of  transitions  between  his  energy  levels  by 
absorbing  or  releasing  quanta  of  the  field.  If  the  energy  of  one  quantum  is  larger  than  the 
typical  level  spacing,  and  provided  that  the  width  of  one  such  one-photon  transition  is 
not  too  large,  then  ,  during  its  time  evolution,  the  system  will  significantly  populate  only 
narrow  resonant  zones  around  the  levels  that  can  be  reached  by  absorbing  or  releasing 
an  integer  number  of  quanta  starting  from  the  initially  excited  levels.  The  motion  will 
thus  resemble  that  of  a  particle  on  a  discrete  lattice,  each  site  in  the  lattice  representing 
a  particular  quasi-resonant  level.  A  jump  from  one  site  to  another  will  correspond 
to  absorption  or  release  of  some  photons.  If  the  levels  were  exactly  resonant  with  one 
another,  tiie  sites  should  be  evenly  spaced  in  energy  by  exactly  one  quantum  and  should 
make  up  a  periodic  lattice,  but  this  will  not  be  the  typical  case.  In  general,  every  quasi 
resonant  level  will  have  a  detuning  from  exact  resonance;  therefore  all  the  sites  in  the 
resonant  lattice  will  be  more  or  less  displaced  from  the  ideal  resonant  configuration  and 
will  make  ip  a  sort  of  a  disordered  lattice. 

This  picture  of  the  dynamics  of  the  dtiven  quantum  systems  is  oversimplified  but 
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it  brings  to  the  light  the  origin  of  dynamical  localization.  Indeed,  provided  that  the 
sequence  of  detunings  can  be  assumed  to  be  irregular  enough,  a  formal  identification 
with  models  of  the  Anderson  type  is  possible.  In  that  case,  no  unlimited  escape  along  the 
quasi-resonant  lattice  is  to  be  expected  -  this  possibility  will  be  ruled  out  by  interference. 

2-  STATEMENT  OP  PROBLEMS 

Is  this  quantum  suppression  a  general  phenomenon? 

At  the  core  of  the  so-called  "dynamical  localization  theory”  lies  the  idea  that  an 
inhibition  of  diffusive  excitation  due  to  quantum  interference  should  be  expected  in 
wide  generality,  upon  quantizing  a  classical  system  in  which  a  chaotic  diffusion  takes 

i 

place.  This  surmise,  supported  by  extensive  numerical  simulations  and  by  some  heuristic 
arguments,  led  to  the  physically  important  prediction,  that  the  quantum  suppression  of 
classical  chaos  should  be  observed  also  in  the  case  of  a  H-atom  in  a  microwave  field.2,3 
This  .vould  result  in  significantly  smaller  ionization  rates  than  classically  expected.  The 
problems  we  addressed  in  our  investigations  were  therefore 

(i)  under  what  conditions  -  if  any  -  can  the  localization  effect  be  observed  in  the  dynamics 
of  real  hydrogen  atoms  in  microwave  fields;  being  aimed  at  decteeting  the  effect  in 
laboratory  experiments,  this  part  of  our  program  involved  a  careful  numerical  simulation 
of  the  actually  realizable  experimental  conditions; 

(ii)  how  can  the  localization  effect  be  circumvented,  giving  rise  to  a  mechanism  of  fast 
unbounded  excitation  eventually  leading  to  ionization, 

(iii)  how  far  the  similarity  between  the  dynamical  localization  and  the  Anderson  local¬ 
ization  can  be  pushed,  and  to  what  extent, can  already  known  theoretical  results  about 
the  latter  be  used  to  describe  the  dynamics  of  classically  chaotic  systems  such  as  atoms 
or  molecules  in  radiation  fields. 


3.-THE  HYDROGEN  ATOM  IN  A  MICROWAVE  FIELD 


The  problem  of  microwave  ionization  of  hydrogen  atoms  has  played  an  extremely  im¬ 
portant  role  in  the  development  of  "quantum  chaology”,  because  it  offers  a  unique 
opportunity  for  experimentally  investigating  the  relevance  of  classical  chaos  in  quantum 

i 

mechanics.  The  behavior  of  the  atom  in  the  microwave  field  can  be  theoretically  and  nu¬ 
merically  analyzed  on  a  very  simple  model  of  a  periodically  driven  nonlinear  oscillator. 
In  the  classical  model  the  onset  of  chaos  provides  a  very  efficient  ionization  mechanism: 
for  any  initial  state  of  the  atom  there  is  a  critical  field  strength  above  which  the  clas¬ 
sical  trajectories  become  chaotic  and  the  atom  starts  absorbing  energy  in  a  diffusive 
way.  This  diffusive  excitation  eventually  leads  to  ionization.  The  theory  of  dynamical 
localization  has  been  applied  to  the  problem  of  a  hydrogen  atom  in  a  microwave  field 
and  has  led  to  an  interesting  prediction,  that  has  been  successfully  checked  in  labora¬ 
tory  experiments,  as  explained  below.2'3,4  Though  in  some  appropriate  frequency  range 
the  process  of  ionization  of  the  atom  is  just  as  strong  as  the  classical  chaotic  picture 
would  predict,  at  higher  frequency  quantum  interference  comes  into  play  and  inhibits 

the  ionization  process.  , 

: 


3.1.-  The  1-d  model  for  microwave  excitation 

The  most  widely  studied  model  for  the  behaviour  of  H-atoms  in  microwave  fields  is  that 
of  an  electron  moving  on  the  half  line  x  >  0  (parallel  to  the  direction  of  the  electric 
field).  The  Hamiltonian  (  in  atomic  units  )  is: 


H  =  — - i-  ex  cos  ut 

2  x 


(i) 


and  is  supplemented  with  a  condition  of  elastic  reflection  at  x  =  0.  In  principle,  the  Id 
model  is  justified  only  for  special  choices  of  the  initial  states  (  extremal  Stark  levels  ), 
but  numerical  simulations  of  the  model  have  given  surprisingly  good  data  for  the  onset 
of  ionization  in  3d  experiments5.  An  equivalent  Hamiltonian  in  a  different  gauge  is 

H  —  g(p  —  <u>-1  sinu <t)2  —  -  (2) 

2  x 
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The  unperturbed  classical  dynamics  (e  =  0)  possesses  a  bound-state  region  in  phase 
space,  and  the  border  of  this  region  is  a  separatrix  (trajectory  of  zero  energy).  Inside 
the  bound-state  region,  action-angle  variables  (n,8)  can  be  introduced.  The  frequency 
fl  of  the  unperturbed  motion  depends  on  the  action  n  according  to  fl(n)  =  7i-J  and 
approaches  0  on  the  separatrix. 

The  application  of  the  resonance  overlap  criterion6  shows  that  for  any  nonzero  field 
strength  e  a  stochastic  layer  exists  close  to  the  separatrix.  Inside  this  layer  the  orbits 
are  chaotic  and  diffuse  until  they  ionize.  The  thickness  of  the  layer  depends  on  e  and  w. 
The  condition  for  strong  ionization  starting  from  a  given  bound  state  is  just  that  the 
stochastic  layer  extends  so  far  downwards  as  to  include  the  given  state.  Given  u>  and 
the  initial  action  n,  the  critical  field  strength  required  for  that  can  be  estimated  an,' 
defines  the  classical  chaotic  threshold 

3.2-The  Kepler  Map 

Following  a  well  known  method  in  the  theory  of  dynamical  systems,  the  continuous 
time  dynamics  ruled  by  the  Hamiltonian  (1,2)  can  be  analyzed  by  means  of  appropriate 
Poincare’  sections,  that  are  introduced  as  follows.  First,  the  time-dependent,  1-freedom 
problem  is  transformed  into  a  time  independent,  2-freedoms  problem  by  introducing  a 
new  couple  of  conjugate  variables  N,<f>  with  <j>  =  cut  (  the  phase  of  the  field  ).  The 
appropriate  Hamiltonian  is  then 

H(p,x,N,<j>)  =  i(p-«u-1sincf)2  -  -  +u>N  (3) 

Z  X 

The  motion  develops  on  a  3d  hypersurface  hi  =  const.  =  A.  From  (3)  it  is  seen  that 
10N  =  A  —  E  where  E  is  the  unperturbed  energy.  The  2d  variety  o  defined  by  p  =  0  can 
now  be  used  to  construct  a  Poincare’  map.  An  orbit  leaving  from  a  point  P  of  o  will 
in  general  cross  again  <r;  the  first  such  intersection  P'  is  taken  as  the  image  of  P  under 
the  map.  Since  N,4>  can  be  used  as  canonical  coordinates  on  <r,  the  map  will  specify 
the  changes  of  N,<f>  between  subsequent  crossings  of  the  variety  a.  If  ew-i  <  1,  these 
crossings  occur  when  i  is  close  to  0,  i.e.,  wher  the  electron  is  close  to  the  aphelion. 
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Therefore,  roughly  speaking  the  map  yields  the  change  in  energy  between  subsequent 
passages  at  the  aphelion,  i.e.  it  corresponds  to  a  stroboscopic  observation  of  the  orbit 
at  intervals  of  time  approximately  equal  to  the  (variable)  orbital  period. 

An  approximate  explicit  expression  of  the  map  can  be  easily  found5’7.  This  ex¬ 
pression  is  especially  simple  if  the  field  frequency  is  much  larger  than  the  unperturbed 
Kepler  frequency  (  i.e.,  if  un*  S>  1  )  *: 

N  =  N  +  k  sin  <j> 

(4) 

4,  =  <j>  +  2™[2(wN  -  A)]"5/2 

Here  k  =  0.822jreai-5/3.  The  above  map  is  defined  only  if  10N  —  A  >  0,  i.e.,  only  for 

negative  energies.  The  iteration  of  the  map  must  be  stopped  when  the  energy  becomes 

positive,  because  then  the  orbit  has  ionized.  Moreover,  the  map  (4)  is  the  result  of  a 

perturbative  calculation  and  is  valid  only  for  not  too  large  k  (see  ref.3). 

Upon  linearizing  the  2nd  eq.(4)  around  some  value  of  N,  corresponding  to  a  value 

Eo  —  — (2«o)_1  of  the  energy,  and  dropping  an  unessential  constant,  one  obtains  the 

Standard  Map:  _ 

N  =  N  +  k  sin  ^ 

(5) 

4>  =  <f>  +  TN 

where  T  =  6jru>2n5.  This  is  a  very  important  result:  for  scaled  frequency  u m3  >  1, 
the  classical  excitation  process  is  locally  described  by  the  Standard  Map.  The  resulting 


picture  of  chaotic  excitation  is  as  follows;  the  electron  mcves  almost  undisturbed  along 
its  unperturbed  orbit  -  apart  from  small  field-induced  oscillations  around  it,  that  do 
not  change  its  average  energy  -  until  it  comes  close  to  the  nucleus.  There  it  receives  a 
’’kick”  described  by  the  Standard  Map.  Above  the  chaotic  threshold,  these  kicks  come 
at  random,  and  the  electron  performs  a  kind  of  a  random  walk  in  energy,  until  it  ionizes. 

3.3-  The  Quantum  Kepler  Map.- 

The  quantization  of  the  Floquet  Hamiltonian  (3)  is  readily  accomplished  by  the  standard 
rule  AT  -»  -ig^.  Once  quantized,  the  variable  N  will  take  only  integer  values;  recalling 

*  In  the  present  form  of  the  map,  the  phase  <f>  15  taken  at  the  perihelion 
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that  —  ioN  coincides  with  the  unperturbed  energy,  apart  from  an  additive  constant,  one 
is  naturally  led  to  identify  N  with  the  number  of  photons. 

In  the  same  vein,  one  can  undertake  a  formal  quantization  of  the  Kepler  Map  (4) 
and  investigate  the  nature  of  the  discrete-time  dynamics  defined  by  the  Quantum  Kepler 
Map  obtained  in  this  way5.  One  finds  that  the  Quantum  Kepler  Map  is  given  by  an 
operator  acting  on  state  vectors  V ’(<j> )  according  to: 

I 

=  e~iil°Pcikm’U(<l>) 

where  , 

*  r  -I — 

Ho  =  2x  |-2u)(fVo  +  A7)j 

the  operator  P  is  the  projector  on  bound  states  N  <  -N0  =  n0/2u0-  This  quantum 

dynamics  can  be  formally  described  as  a  motion  on  a  discrete  lattice  where  the  sites  are 

labelled  by  the  integer  values  of  N.  This  has  the  following  physical  meaning.  Starting 

i 

from  some  initial  hydrogenic  bound  state,  and  by  absorbing  or  releasing  photons,  the 
atom  will  jump  from  one  state  to  another  in  a  ladder  of  ’quasi  resonant’  bound  states 
("photonic”  states).  These  states  correspond  to  the  sites  in  the  lattice,  and  the  electron 
can  jump  from  one  site  to  another  by  absorbing  or  emitting  some  photons. 

Since  the  classical  map  describes  the  motion  only  in  the  bound  state  region,  from 
which  the  electron  may  leave  when  ionization  occurs,  tl:e  quantum  map  can  account 
only  for  the  absorption  of  a  maximum  number  of  photons,  given  by  the  number  of 
photons  required  to  ionize:  Ni  =  (2nj<v)  .  In  other  words,  the  quantum  map  defines 

the  dynamics  on  a  finite  lattice.  The  quantum  map  also  includes  a  dissipation  which 
describes  direct  ionization:  this  is  active  only  on  a  few  photonic  states  close  to  the 
1-photon  ionization  border. 

As  can  be  guessed  from  the  similarity  of  the  classical  Kepler  map  to  the  Standard 
Map  (5),  the  Quantum  Kepler  Map  is  closely  related  to  the  Quantum  Kicked  Rotator, 
that  will  be  presently  quickly  reviewed. 
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3.5  -Localization  in  the  Quantum  Kicked  Rotator 


The  quantum  kicked  rotator'is  obtained  by  quantizing  the  Standard  Map  (5).  Its  motion 
can  be  formally  described  as  a  quantum  motion  on  the  discrete,  infinite  lattice  defined 
by  the  integer  values  of  N  (  that  has  now  the  meaning  of  angular  momentum  ).  This 
motion  is  always  localized:  a  wave  packet  initially  concentrated  at  some  site  No  will 
start  spreading  on  the  lattice,  but  after  a  while  it  will  "freeze”  into  an  average  steady- 
state  configuration.  The  corresponding  probability  distribution  over  the  sites  will  be 
fluctuating  around  an  average  distribution: 


m  = 


1+ 


2\N  —  iV0 1 

l 


exp  - 


2|JV  -  jVo| 


l 


(6) 


where  l  is  the  localization  length.  In  the  semiclassical  regime,  and  in  the  region  of  strong 
classical  chaos,  one  has  the  remarkable  estimate  l  ~  |fc2. 

The  above  estimate  comes  from  a  very  general  argument  that  can  be  summarized  as 
follows11.  Let  us  consider  a  periodically  perturbed  quantum  system  which  has  a  chaotic 
classical  limit.  Given  the  states  ifi(k)  of  the  system  after  k  =  1,2,..,,  A1  periods  of  the 
perturbation,  we  can  investigate  the  quasi-energy  spectrum  by  performing  a  discrete- 
time  Fourier  analysis  of  ther  string  c*  =  (ij’(0)\tj>(k))  where  (.|.)  is  the  scalar  product. 
Given  Cfc  up  to  k  =  N,  we  can  resolve  the  spectrum  on  a  scale  ~  jj.  On  the  other 
hand  we  can  expect  the  number  of  frequencies  significantly  contributing  in  the  motion 
up  to  time  N  to  be  on  the  same  order  as  (he  number  of  unperturbed  levels  excited  up 
to  the  same  time;  this  number  is  roughly  An  =s  (DN)l  with  D  the  classical  diffusion 
coefficient.  Therefore,  on  the  scale  1/N  tile  quasi  energy  spectrum  will  be  resolved  by 
just  (RiV)>  values.  This  means  that  for  sufficiently  large  N  gaps  will  appear  in  the 
spectrum,  roughly  for  N  s:  (DN)i,  i.e.,  for  N  &  D.  For  such  ~N  recurrences  should  be 
expected  to  appear  and  to  slow  down  the  excitation.  The  spread  An  of  the  wave  packet 
at  the  same  lime  N  will  be  on  the  order  of  *  he  maximum  spread  attainable  by  the  wave 
packet  in  the  course  of  its  evolution;  since  the  latter  is  just  given  by  the  localization 
length,  we  get  l  =3  D  ss  k2/ 2  where  the  latter  estimate  is  well  known  and  comes  from 
estimating  the  classical  diffusion  coefficien.  under  a  random  phase  assumption  that  is 
very  well  justified  in  the  case  of  strong  classical  chaos. 


tir 
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3.5  -Localization  in  the  H-atom 


The  dynamical  localization  theory  for  the  Hydrogen  atom  in  a  microwave  field  is  basi¬ 
cally  an  extrapolation  to  that  problem  of  the  resells  obtained  for  the  qnantnm  rotator 
(  QKR  ),  a  bridge  bring  provided  by  the  Qnantnm  Kepler  Map  (  QKM  ).  Tins  extrap¬ 
olation  rests  on  some  assumptions: 

-  that  the  QKM  is  in  some  sense  an  approximation  for  the  real  Schrodinger  equation. 
This  is  by  no  means  obvions:  the  connection  between  the  discrete  time  defined  by  the 
number  of  iterations  of  the  Map  and  the  real  time  is  quantum-  mechanically  lost  as 
soon  as  the  wave  packet  spreads  over  states  which  have  significantly  different  classical 
orbital  periods. 

-  that  the  similarity  in  structure  between  the  QKR  and  the  QKM  justifies  the  application 
to  the  QKM  of  known  results  about  the  OKR.  The  essential  difference  between  the 
QKM  and  the  QKR  is  that  the  former  lives  in  a  finite  lattice,  and  that  this  lattice 
has  a  ’leaky  end’  due  to  direct  ionization  from  photonic  states  close  to  the  ionization 
border.  Nevertheless,  as  long  as  the  wave  packet  does  not  significantly  populate  these 
border  states,  the  dissipation  can  be  neglected  and  the  QKM  and  QKR  dynamics  will 
be  similar  -  in  their  gross  features  at  least.  This  will  certainly  be  the  case  if  l  <C  A’r, 
because  the  wave  packet  will  "freeze”  before  haring  a  chance  to  significantly  populate 
the  border  sites.  Since  as  explained  above  l « Jfc2/2  for  the  KR  and  the  "kick  strength” 
of  the  Kepler  Map  is  k  =  0.0822wew-5/3,  one  should  expect  the  localization  length  in 
the  number  ol  photons  to  be  given  by 

l  =  3.33eJa>-10/5 


which  is  well  confirmed  not  only  by  the  numerical  simulation  of  the  quantum  Kepler 
Map  but  even  by  the  numerical  solution  of  the  Schrodinger  equation  (Figs.l),  thus  fully 
supporting  the  validity  of  the  Kepler  Map  approach  also  for  the  quantum  problem. 
Given  this  explicit  formula  for  /,  the  condition  for  localization  /  <  Nj  takes  the  form: 

*/« 


e  <£-  = 


(6.6«0)1/2 


(7) 
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There  is  a  wide  range  of  parameters  where  ibis  condition  can  be  satisfied  while 
keeping  the  field  strength  well  above  the  classical  chaos  border.  In  this  region,  a  phe¬ 
nomenon  of  dynamical  localization  quite  similar  to  the  one  observed  in  the  QKR  will 
prevent  the  atom  from  ionizing  so  fast  as  if  should  classically.  The  ”  frozen”  distribution 
will  have  the  form  (6)  with  N  the  number  of  absorbed  photons. 

3.5-Delocalization 

If  l  is  increased  (  c-g.,  by  increasing  the  field  strength  )  while  keeping  Nj  constant,  the 
above  picture  will  lose  its  validity.  Qualitatively  one  can  expect  the  ionization  rate  to 
increase  due  to  the  onset  of  a  diffusive  excitation  more  or  less  like  the  classical  one. 
Conditions  for  this  are  that  l  Nj  and  that  at  the  same  time  k  C  Nj  (  implying  that 
a  large  number  of  kicks  are  needed  to  ionize  ). 

3.8-  Numerical  Checks 

The  localization  phenomenon  for  the  H-atom  in  a  microwave  field  has  been  confirmed 
by  extensive  numerical  simulations  of  the  time-dependent  Schrodinger  equation  since 
19842.  In  addition  to  that  the  Kepler  map  yields  an  indication  of  great  heuristic  value, 
namely  that  the  shape  of  the  localized  distribution  should  be  exponential  in  the  num¬ 
ber  of  absorbed  photons  (  formula  (6)  ).  This  has  also  been  confirmed  by  numerical 
computations3.  In  fig  1  an  example  is  shown  of  localized  distributions  obtained  from 
the  numerical  integration  of  the  Schrodinger  equation.  The  distribution  over  the  un¬ 
perturbed  hydrogenic  states  (labelled  by  the  quantum  number  n),  ,  suitably  averaged 
in  time,  is  plotted  (  in  logarithmic  scale  )  versus  the  number  of  absorbed  photons 
N  =  (2n2)-1  —  (2nJ)_1  where  no  =  100,200  is  the  quantum  number  of  the  initially 
excited  state.  Here  tonj  =  3.5  .The  exponential  character  of  the  distribution  is  fairly 
evident,  as  well  as  the  role  of  the  quasi  resonant  ’photonic’  states. 
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3.0-  Beyond  the  Id  Model 

Although  the  above  theory  was  developed  for  the  1-dimensional  model,  its  predictions 
turn  out  to  be  valid  even  for  more  realistic  situations.  As  we  already  mentioned,  the 
experimentally  observed  ionization  thresholds  could  be  satisfactorily  reproduced  by  nu¬ 
merical  simulation  of  Id  classical  models;  on  qualitative  grounds,  it  was  held  that  almost 
one  dimensional  orbits,  i.e.  orbits  very  elongated  in  the  field  direction  could  be  more 
easily  ionized  so  that  they  would  give  a  dominant  contribution  to  the  ionization  rate  in 
the  threshold  region. 

We  have  investigated  a  two-dimensional  model  in  which  the  magnetic  quantum 
number  with  respect  to  the  field  direction  is  0.  We  were  able  to  find  a  Kepler  Map 
that  approximately  describes  the  changes  of  the  canonical  variables  between 

subsequent  passages  at  the  aphelion;  C  is  angular  momentum,  and  ip  its  conjugated 
phase.  Although  this  map  is  quite  complicated3  ,  an  approximated  analysis  shows  that 
the  evolution  of  elongated  orbits  is  essentially  one-dimensional  -  i.e.,  it  develops  as 
predicted  by  the  Id  model  -  due  to  the  existence  of  an  approximate  integral  of  motion. 
Over  not  too  long  times  (  on  the  order  of  the  experimental  interaction  times  )  the 
two  degrees  of  freedom  behave  as  if  they  were  uncoupled.  This  theoretical  argument  in 
support  of  the  validity  of  the  Id  model  even  in  2d  situations  was  confirmed  by  numerical 
simulations  of  a  quantum  2d  atom,  that  yielded  evidence  of  an  exponential  localization 
in  the  number  of  photons  quantitatively  similar  to  the  one  occurring  in  the  Id  model 
(Fig-2). 


3.8-Experiniental  Verifications.- 

Experimental  results  on  microwave  ionization  of  highly  excited  Hydrogen  atoms  in  the 
region  wnj  >  1  were  reported  for  the  first  time  in  1988  by  Galvez  et  at3 .  Until  then,  nu¬ 
merical  computations  and  experiments  in  .lie  region  umj  <  1  had  shown  that  classical 
mechanics  accounts  very  well  for  most  of  tre  experimental  results.  The  1988  exper¬ 
imental  data  clearly  show  that  in  the  high-frequency  region  the  quantum  ionization 
thresholds  rise  above  the  classical  ones.  These  data  provide  experimental  evidence  for 
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the  quantum  suppression  phenomenon.  This  was  confirmed  by  new  experiments  by  an¬ 
other  experimental  group  in.  1989;  in  this  case,  the  experimental  results  were  also  found 
to  agree  with  some  quantitative  predictions  of  localization  theory10. 

A  typical  microwave  ionization  experiment  produces  "ionization  curves"  like  those 
shown  in  Fig.3,  which  displays  some  of  the  experimental  data  by  Galvez  et  al.  For  given 
oj  and  n0  the  field  strength  is  increased  until  lOpercent  ionization  is  achieved.  The 
corresponding  value  of  the  ’scaled  field’  Co  =  is  plotted  versus  the  ’scaled  frequency’ 
u0  =  uni;  the  latter  is  varied  by  changing  n0  .  The  data  in  Fig.3  were  obtained  with 
u  =  36.02  GHz.  Fig.3  shows  that  for  u0  >  2  the  lOpercent  experimental  thresholds  ( 
circles  )  lie  systematically  above  the  (  numerically  computed  )  classical  thresholds  (  the 
crosses  in  Fig.3  ). 

A  crucial  remark  is  that  the  experimentally  measured  ionization  does  not  only  include 
those  electrons  that  have  escaped  into  the  continuum,  but  also  those  excited  above 
some  cutoff  level  nc  (  for  a  detailed  discussion  see  ref.5).  Thanks  to  this  fact,  the 
localization  theory  can  sometimes  produce  a  quantitative  estimate  for  the  thresholds. 
This  is  the  case  when  only  a  small  part  of  the  total  probability  flows  into  continuum 
in  the  given  interaction  time.  In  that  case,  (i)  the  ’ionization  probability’  P  will  be 
essentially  given  by  the  total  probability  on  the  bound  states  above  the  cutoff  nc,  and 
(ii)  the  final  distribution  over  the  photonic  states  will  be  of  the  form  (6).  Therefore,  P 
can  be  computed  from  the  distribution  (6).  In  principle,  this  procedure  can  be  reversed 
in  order  to  find  what  the  field  strength  must  be,  in  order  that  P  —  0.1.  This  yields 
a  simple  formula  11  which  corresponds  to  the  smooth  line  in  Fig.3.  More  precisely,  if 
we  define  the  ionization  probability  as  the  total  probability  above  a  fixed  level  n,  the 
threshold  for  10%  ionization  is  obtained  from  the  condition 


0.01 


/•: 


f(N)dN 


witli  N  =  (rij  2  —n~2)/2u.  The  rhs  of  the  last  eqn.  can  be  computed  as  a  function  of 
i  by  using  eqn(6),  and  the  resulting  equation  can  be  solved  for  e. 

However,  in  order  that  this  procedure  be  a  self-consistent  one,  the  field  strength 
and  the  localization  length  found  in  this  way  must  not  be  so  large  as  to  significantly 
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populate  tlie  fast-ionizing  states  dose  to  the  ionization  border;  otherwise,  the  dissipation 
cannot  be  assumed  to  be  small.  In  such  a  case  the  ’localization  formula’,  which  does 
not  account  for  continuum,  cannot  be  used  -  it  would  be  found  to  overestimate  the 
actual  thresholds.  For  instance,  this  will  happen  if  the  cutoff  nc  is  pushed  very  far  from 
the  initially  excited  state  -  it  will  not  be  possible  in  that  case  to  transfer  a  lOpercent 
probability  on  the  bound  states  beyond  the  cutoff,  without  populating  the  border  states. 
One  such  case  that  lies  beyond  the  limits  of  applicability  of  the  localization  formula  was 
reported  in  refs  5,9. 

A  further  confirmation  of  the  localization  effect  was  provided  by  experimental  work 
of  J.Bayfield  and  D.W.Sokol10.  This  work  was  especially  designed  in  order  to  check  the 
predictions  of  localization  theory  and  was  carried  out  in  close  parallel  with  numerical 
simulations  of  the  experimental  setup.  Wheras  in  the  case  of  Galvez  ei  al.  atoms 
were  prepared  in  a  microcanonical  distribution  of  states  with  a  given  principal  quantum 
number,  in  Bayfield’s  case  the  magnetic  quantum  number  along  the  direction  of  the 
field  could  be  assumed  to  be  zero;  the  situation  here  was  therefore  ’’two-dimensional” 
as  compared  to  the  ’’three-dimensional”  situation  of  Galvez  et  al.  Another  point  of 
difference  was  that  the  microwave  interaction  took  place  inside  a  wave  guide  (  Galvez  et 
al  had  a  microwave  cavity  instead  )  where  atoms  were  exposed  to  a  microwave  pulse  of 
frequency  varying  in  a  band  12.4  to  18  GHz,  with  a  nominal  envelope  sin(t/7.!>nsec). 
Finally,  a  static  field  0.87  Volt/cm  was  superimposed  to  the  microwave  field;  this  field 
would  have  produced  static-field  ionization  of  states  with  principal  quantum  number 
164,  considerably  higher  than  the  adjustable  cutoff  level  n  used  to  empirically  define 
ionization. 

Figure  (4)  shows  thresholds  for  10%  ionizations  obtained  experimentally,  by  quan¬ 
tum  numerical  simulations  on  the  Id  model,  and  by  the  above  sketched  theoretical 
formula.  (  solid  line  ).  To  be  emphasized  is  that  the  theory  does  not  describe  either 
quantum  or  classical  resonance  effects,  as  it  is  based  upon  only  the  smoothed  behaviour 
of  stationary  state  distributions  in  quantum  number. 
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4.-  HIGHER  DIMENSIONAL  LOCALIZATION 


How  can  the  localization  effect  be  circumvented?  In  the  theory  of  Anderson  localization 
in  disordered  solids  this  question  admits  of  a  variety  of  different  answers,  depending 
on  the  spatial  dimension  of  the  solid  themselves.  In  the  one  dimensional  case,  i.e.  for 
a  solid  shaped  like  a  bar  or  a  wire,  the  localization  becomes  ineffective  sis  soon  sis  the 
localization  length  becomes  larger  than  the  length  of  the  solid.  The  delocalization  in 
the  H-atom  problem  has  the  same  origin:  here  the  length  of  the  ’’solid”  is  defined  by  the 
number  of  photons  that  will  carry  the  electron  into  the  continuum  (  see  the  discussion 
in  Sec.1.2)  ,  and  the  condition  that  this  length  be  smaller  than  the  localization  length 
is  just  eqn.(7).  For  two-  and  three-dimensional  solids,  besides  this  obvious  finite-sample 
effect,  other  mechanisms  exist  that  make  the  the  effect  of  Anderson  localization  milder. 
It  is  a  natural  question  whether  these  higher-dimensional  localization  effects  have  a 
dynamical  counterpart,  too. 

As  summarized  above,  the  dynamics  of  a  quantum  system  subjected  to  a  perturba¬ 
tion  periodic  in  time  is  formally  ’’conjugate”  to  a  one-dimensional  localization  problem. 
It  turns  out  that  higher  dimensional  localization  problems  have  a  dynamical  counter¬ 
part  in  the  dynamics  of  systems  subjected, to  a  number  of  incommensurate  frequencies. 
Notice  that  the  Kicked  Rotator  is  subject  to  infinitely  many  frequencies  commensurate 
frequencies  given  by  the  integer  multiples  of  the  kick  frequency  and  that  the  hydrogen 
atom  is  subject  to  a  monochromatic  field.  It  is  known[4j  that  the  addition  of  a  second 
incommensurate  frequency  in  the  kicked-rotator  model  does  not  destroy  the  localization 
effect;  nevertheless,  the  localization  lengths  are  sharply  increased  in  comparison  to  the 
one  frequency  case,  in  a  way  that  somehow  reproduces  known  results  from  the  tliory  of 
2-d  disordered  solids  (’’films”). 

The  most  interesting  case  is  the  one  with  three  incommensurate  frequencies.  As 
we  shall  presently  explain,  this  corresponds  to  a  three  dimensional  localization  prob¬ 
lem.  According  to  the  theory  of  Anderson  localization  a  transition  from  localized  to 
extended  states  should  be  expected  as  the  coupling  strength  is  varied.  The  reason  why 
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investigating  the  possibility  of  a  similar  transition  in  a  dynamical  model  is  important  is 
twofold.  First,  one  needs  to  ascertain  whether  the  localization  effects  that  prevents  dif¬ 
fusive  excitation  can  be  drcumventerd,  so  that  in  some  cases  at  least  quantum  systems 
can  prove  as  unstable  as  their  classical  counterparts  are.  Second,  the  precise  details 
of  the  Anderson  transition  are  not  yet  completely  understood;  the  theoretical  analysis 
is  complicated  and  the  numerical  simulation  prohibitive.  In  contrast,  the  numerical 
analysis  of  dynamical  models  is  much  less  expensive;  as  a  matter  of  fact,  in  this  way 
we  were  able  to  obtain  some  interesting  indications  at  a  relatively  low  computational 
cost12. 

4.1-  Anderson  Transition  in  a  Dynamical  Model 

We  studied  the  motion  of  a  particle  on  a  circle,  described  by  the  time-dependent  Hamil¬ 
tonian: 

+CO 

H  =  H0+V(6,t)  Y,  «(<-*)  (8) 

S—~  CO 

The  2nd  term  describes  kicks  occurring  periodically  in  time  with  period  one.  The  free 
evolution  between  kicks  is  given  by  the  Hamiltonian  Ho'- 

Ho\n)  =  En\n),  |n)  =  (2  n)'1'2  exp(in9) 

We  assume  the  eigenvalues  En  to  be  random  numbers  uniformly  distributed  in  (0,2rr). 
We  also  assume  V  to  explicitly  depend  on  time  according  to 

V  =  V(Q,9\  +uJit,^2  +  u>2t) 

with  V  a  periodic  function  of  its  three  arguments  to  be  specified  later,  and  0,  l2  arbitrarily 
prescribed  phases.  Since  we  would  like  u>,,2  to  be  incommensurate  with  each  other 
and  also  with  the  frequency  of  the  kicks,  we  take  <vj  =  2ttA_1  and  u>2  =  2xA~2  with 
A  =  1.3247...  the  real  root  of  the  cubic  equation  xz  —  x  —  1  =  0.  With  such  a  choice, 
uq,2  are  a  ’most  incommensurate”  pair  of  numbers.  The  evolution  of  the  system  from 
just  after  one  kick  to  just  after  the  next  is  given  by  (  we  take  here  ft  =  1  ): 

+  1)  =  e-1  “'•<+,>e-H°V’((M)  (9) 
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Expressed  in  this  form,  the  rotator  dynamics  is  very  easily  simulated  because  the  time 
dependence  of  V  is  explicitly  known. 

We  shall  now  elucidate  how  is  this  model  connected  with  a  3d  tight-binding  model 
of  the  Anderson  type.  First  of  all  we  consider  the  phases  6j  l2  as  new  dynamical  variables, 
with  conjugate  momenta  nli2.  Then  we  consider  the  Hamiltonian 

+  OO 

H1  =  Ho{h)  +  +  wziii  +  'y  ^  —  s)  (10) 

<=—00 

with  rii|2  =  —id/dd i,j.  The  last  hamiltonian  describes  a  rotator  with  three  freedoms 
0,61,62  subjected  to  periodic  kicks,  the  strength  of  which  is  no  more  time  dependent. 
The  one-period  propagator  for  this  rotator  is  now 

£«-iV(0,St,ea)g— t(Ho(n)+w,ni+W2A,)  (ll) 

In  order  to  recognize  that  the  3d  quantum  model  defined  by  (10)  and  the  Id  model 
defined  by  (8)  are  substantially  equivalent,  it  suffices  to  rewrite  the  Schrodinger  equation 
for  the  3d  model, 

ijl'K6,8u62,i)  =  H,W,eu62,t) 
in  the  interaction  representation  defined  by 
^(6,6i,02,t)  = 

In  this  way,  one  immediately  recovers  the  Schrodinger  equation  with  Hamiltonian  (8) 
for  the  wave  function  ip. 

We  can  now  apply  to  the  3-rotator  (10)  a  well-known  transformation1  that  will 
map  it  onto  a  3d  tight-binding  model.  Thanks  to  this  transformation,  the  problem 
of  determining  the  quasi-energy  eigenvalues  and  eigenvectors  of  the  3-rotator  will  be 
turned  into  the  problem  of  solving  the  equation: 

TnUu  *1-  y  ^  Wrun+r  =  tu„  (12) 

where  n  =  (n, ,  n2)  and  r  label  sites  in  a  3d  lattice, 

T„  =  -tan[0.5(E„  f  n,u,  -f  n2iv2  +  A)) 
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Here  X  is  quasi-energy,  Wr  are  the  coefficients  of  a  threefold  Fourier  expansion  of 
tan[0.5V(6,6i,62)]  and  e  =  —Wq.  We  now, choose 

V(8, 61,62)  —  —2arctan[2k(cos8  +  cos  (?j  +  cos 62) 
and  eqn.(12)  becomes  1 

Tuun  +  k  ^2  ut  =  0 

r 

where  the  sum  includes  only  nearest  neighbours  to  n. 

After  all  these  formal  manipulations,  the  original  rotator  problem  has  been  turned 
into  a  tight-binding  model  of  the  Anderson  type;  the  nature  of  the  eigenstates  of  the 
latter  will  determine  the  nature  of  the  dynamics  of  the  rotator,  that  will  be  localized  or 
delocalized  (  in  the  sense  of  unbounded  excitation  taking  place  )  according  to  whether 
these  eigenfunctions  are  localized  or  extended.  The  disorder  in  the  model  (12)  is  given 
by  the  pseudorandom  character  of  the  potential  T„. 

The  model  was  investigated  by  numerically  simulating  the  dynamics  of  the  rotator, 

I 

starting  with  the  same  initial  state  and  with  increasing  values  of  the  coupling  constant 
k.  A  transition  between  two  types  of  motion  was  detected  around  a  value  kcr  ~  0.47, 
the  dynamics  being  localized  for  k  <  kcr  and  delocalized  for  k  >  kcr.  In  the  former  case 
the  (  time-averaged)  steady-state  distributions  were  found  to  be  exponential  (Fig.5)  and 
the  corresponding  localization  lengths  could  be  determined.  In  the  latter,  unbounded 
diffusive  excitation  occurs,  with  Gaussian-like  distribution  of  width  linearly  increasing 
with  time  (Fig.O)  ,  so  that  the  diffusion  coefficients  could  be  determined.  A  remarkable 
result  was  obtained  from  the  analysis  of  the  behaviour  of  the  inverse  localization  length 
7  and  of  the  diffusion  coefficient  close  to  the  transition  point  (Fig.7):  we  found  that 

D  ~  const. |fc  -  bcr|*,  7  ~  const.(\k  -  kcr |"’ 

with  s  -  1.25  and  u  ~~  1.5,  which  are  close  to  values  theoretically  obtained  in  localization 
theory  from  diagrammatic  computations. 

This  result  of  ours  enforces  a  striking  similarity  between  the  dynamics  of  quantum 
systems  under  almost  periodic  perturbations  and  tight-binding  models  with  spatial  di¬ 
mension  given  by  the  number  of  incommensurate  frequencies. 


3.-  CONCLUSIONS  AND.  RECOMMENDATIONS 


The  core  of  the  research  activity  summarized  in  this  final  report  was  the  development 
of  the  dynamical  localization  theory  and  its  application  to  a  specific  problem  in  atomic 
physics.  Originally  motivated  by  a  theoretical  interest  in  the  possible  survived  of  effects 
of  classical  stochasticity  in  quantum  mechanics,  this  research  has  brought  into  the  light  a 
class  of  generic  behaviours  of  small  quantum  systems  subjected  to  periodic  perturbations 
such  as,  e.g.,  atoms  and  molecules  in  external  fields.  The  excitation  of  such  systems  is 
strongly  limited  by  a  coherence  effect  that  is  basically  the  same  that  is  responsible  for  the 
Anderson  localization  in  solid  state  physics;  by  generalizing  some  typical  concepts  and 
methods  of  the  theory  of  Anderson  localization,  we  were  able  to  get  a  simple  albeit  rough 
quantitative  description  of  the  excitation  of  atoms  in  microwave  fields  that  has  stood 
the  test  of  laboratory  experiments.  The  dynamical  localization  theory  has  predicted 
the  quantum  suppression  of  diffusive  ionization  that  was  first  confirmed  by  numerical 
simulations  and  then  observed  in  laboratory. 

Notwithstanding  this  undoubted  success,  much  work  has  still  to  be  done.  Although 
already  in  its  present  form  the  theory  represents  a  statistical,  non-perturbative  approach 
that  can  prove  useful  in  a  variety  of  time-dependent  problems,  more  theoretical  work 
is  needed  in  order  to  completely  understand  how  does  localization  stem  out  of  the 
Schrodinger  equation.  Even  more  important,  the  firmly  established  similarity  of  the 
dynamical  localization  with  the  better  known  Anderson  localization  opens  the  way  to 
further  exciting  developments. 

As  we  have  emphasized  several  times  in  this  report,  the  predictions  of  our  theory  are 
essentially  statistical,  as  a  result,  the  agreement  with  experimental  data  holds  only  ”in 
the  average”.  Indeed  both  the  numerical  and  the  experimental  data  exhibit  more  or  less 
relevant  deviations  from  the  average  theoietical  prediction.  This  is  hardly  surprising, 
because  the  dynamical  localization  theory  was  just  meant  to  yield  a  gross  description  of 
the  quantum  dynamics;  actually,  as  it  was  pointed  out  in  ref3,  one  would  expect  even 
stronger  fluctuations  than  were  actually  observed.  Indeed,  according  to  the  theory, the 
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quantum  distribution  is  strongly  fluctuating  around  its  average  exponential  shape,  and 
these  fluctuations  will  affect, also  the  ionization  rate  which  will  have  a  very  irregular  fine 
structure  . 

Short  of  representing  a  failure  of  our  theoretical  approach,  these  "fluctuations”  are 
an  essential  part  of  the  theoretical  picture  of  localization  theory  and  in  the  case  of  An¬ 
derson  localization  their  statistical  description  is  a  major  goal  of  the  theory  .  They  are 
not  due  to  uncontrolled  changes  in  the  environment  or  to  any  form  of  external  noise; 
they  are  perfectly  reproducible,  so  that  in  principle  the  very  use  of  the  word  "fluctu¬ 
ations”  might  be  questioned.  Though  thlse  irregularities  stem  from  purely  quantum 
coherence  effects  -  nothing  comparable  is  to  be  found  in  the  classical  ionization  rates  -, 
it  may  well  be  the  case  that  the  seemingly  irregular  structures  that  characterize  the  ion¬ 
ization  curves  correspond  to  identifiable  structures  in  the  classical  phase  space  (  "scars” 
).  Nevertheless,  an  individual  description  of  such  local  structures  appears  hardly  feasi¬ 
ble.  Although  there  have  been  isolated  successes  in  associating  some  quantum  details 
with  classical  objects  such  as  periodic  orbits  or  "homoclinic  tangles”,  ,  this  kind  of 
classical  structures  are  of  course  individually  unpredictable,  as  is  for  the  time  being  the 
mechanism  why  some  of  them  find  sometimes  their  way  to  "scarred”  quantum  states. 

Therefore,  from  a  physical  standpoint,  a  statistical  approach  is  recommended, 
aimed  at  quantitatively  describing  how  do  these  irregular  variations  of  ionization  rates 
scale  with  the  relevant  parameters. 

On  more  general  grounds,  these  "fluctuations”  are  but  one  more  instance  of  a  sensi¬ 
tive  dependence  of  quantum  reaction  rates  on  various  parameters  that  has  been  observed 
in  several  areas,  from  nuclear  physics  to  solid  slate.  Some  universal  characteristics  are 
now  emerging,  that  are  conveniently  modelled  by  Random  Matrix  Theory,  and  many 
investigators  of  these  different  fields  are  now  adopting  as  a  working  hypothesis  that  a 
unified  approach  at  such  different,  purely  quantum  "fluctuation”  phenomena  may  be 
based  on  the  chaotic  nature  of  the  underlying  classical  mechanics.  Insofar  as  future 
work  is  concerned,  we  believe  that  just  this  one  is  the  most  promising  indication  of  the 
work  we  have  been  summarizing  in  this  report. 
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5.-  FIGURE  CAPTIONS 

Fig.l-  Probability  distribution  over  the  unperturbed  levels  averaged  bom  450  to  500 
microwave  periods  (full  curve)  versus  the>  number  of  photons  N^,  for  Co  =  0.03, «o  = 
3.5.The  straight  line  results  from  a  least  squares  fitting  of  the  maxima  of  the  distribution 
in  each  photon  interval.  The  dots  give  the  steady  state  distribution  obtained  by  iterating 
the  quantum  map.(a)  no  =  100,  (b)no  =  200.  For  comparison,,  the  classical  distribution 
over  unperturbed  levels  is  also  shown  in  (b)  (dashed  curve)(  From  ref.(3)). 

Fig.2-Classical  (dashed  curve)  and  quantum  (full  curve)  distribution  functions  versus 
the  number  of  absorbed  photons  and  for  the  2-dimensional  model.  The  straight  dotted 
line  is  the  one-dimensional  quantum  theoretical  exponential  distribution;  the  dotted- 
dashed  line  is  the  analytical  solution  of  the  diffusion  equation. 

Fig.3-Scaled  10%  threshold  fields  from  ref.(9).  The  dashed  line  is  the  quantum  theoret¬ 
ical  prediction  according  to  dynamical  localization  theory. 

Fig.4-  A  point  by  point  comparison  of  experimental  and  quantum  mechanical  values  for 
the  microwave  field  strength  for  10%  ionization  probability,  as  a  function  of  microwave 
frequency.  The  field  and  frequency  are  scaled  according  to  ui0  =  unl,F0  =  enj.  The 
theoretical  points  are  shown  as  solid  triangles.  The  experimental  points  arc  connected 
by  a  dashed  curve  drawn  through  the  entire  data  set.  Values  of  the  initial  level  no 
and  of  the  cutoff  level  n  are:  full  circles, 64/1 14;  crosses, 68/114;  open  circles,  71/114; 
squares, 80/120;  triangles,  86/130;  plusses,,  94/130;  diamonds,  98/130.  The  dotted  line 
is  the  classical  chaos  border;  the  solid  curve,  the  quantum  border  . 

Fig.5- Example  of  a  localized  steady-state  distribution  over  the  unperturbed  levels  of  the 
three-frequencies  model,  averaged  over  5000  iterations  m  the  interval  95000(tjl00000. 
Here  k  =  0.38. 

Fig. 6.  Example  of  a  Gaussian  distnbn  ion  over  the  unperturbed  levels  in  the  delocalized 
regime  of  the  three-frequencies  model.  The  distribution  is  averaged  over  1000  iterations 
in  the  interval  19000jt;20000. 
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Fig.7.-  Diffusion  rate  D  (dots)  and  inverse  localization  length  7  as  functions  of  the  cou¬ 
pling  strength  k.  Error  bars  were  obtained  from  statistics  over  ten  different  realizations 
of  the  random  spectrum.  The  dotted  lines  result  from  a  3-parameters  least  squares  fit. 
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